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Abstract 

We study nonlinear effects in the QED ladder Schwinger-Dyson(SD) equation. Without further 
QQ ' approximations, we show that all nonlinear effects in the ladder SD equation can be included in 

^S] , the effective couplings and how a linear approximation works well. The analyses is generalized in 

the case of the improved ladder calculation with the Higashijima-Miransky approximation. 
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The dynamical symmetry breaking (DSB) is widely studied by the Schwinger-Dyson (SD) 
equation. (See references m and H, and the references therein. ) In QED, the lowest-order 
\^ ' approximation is a ladder approximation. It was pointed out that there is DSB of the chiral 

^^ . symmetry in the QED ladder SD equation with a finite cut off. (See, e.g., p, 0, 0]. ) In QCD, the 

0^ ' calculation was modified by introducing a running coupling. The modified calculation is called an 

improved ladder approximation and gives a good results even quantitatively. 
,S^ ' Although the original SD equation is an integral equation, it is sometimes rewritten in the 

Mh, form of a nonlinear differential equation. (See, e.g., 0, 0, nl 0. ) Even in the most simple ladder 

O ' approximation, the nonlinear equation can not be solved exactly and a linear approximation is 

D , widely used. (See, e.g., [§|, ||. ) The aim of this paper is to re-examine the nonlinear effects 

1-^ ' in the QED ladder SD equation and show how the linear approximation works well. After that, 

we generalize our analyses in the case of the improved ladder calculation with the Higashijima- 
Miransky approximation. [[7|, ^ 

In Landau gauge, the QED ladder Schwinger-Dyson equation is given by 



X 



Bip')^3e'J 



Too, (1) 



{2tt)U^ + B^q^) {p - qy 



where B{p^) is a fcrmion self-energy with four momentum p and mo is a bare mass of the fermion. 
We restrict our discussions on Euclidean region only. Putting p = me^{m is a some energy scale) 
and rescaling x{t) = B{p^)/p, after some algebra, we get a nonlinear second differential equation 
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x + 2x+-x + \- ^ = (2) 

A 1 + x^ 



with two boundary conditions 



and 



3t I 

vi = e~{x + -x) -^0 {t^ -oo) (3) 

V3 = e^{x + -x)^ -^ (i->oo). (4) 
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In these equations, A — -^j^, i = j§ and x = ^^f . Since we are interested in DSB of chiral 
symmetry, we put toq = 0. The uhra-violet boundary condition (UVBC) (0) becomes 

o 

V3 = e2(x+-x)^0 (i— >oo). (5) 

The equation (y) can be rewritten as 

dvi xe~ dvs xe^ 

-—- = —A ^ or — — = —A 7T- (6) 

dt 1 + x^ dt l + x^ ^ ' 

We can construct a potential V{x) as [g[ 

xr/ X 3 o A , , ,, dVix) Z ^ X ,„, 

V{x)^-x^ + -\og{l + x^)- -^= x + \-—^ 7 

8 2 da; 4 1 + x^ 

and a hamihonian H{y,x) 



Hiy,x) = ^y' + Vix), (8) 

where y = i; is a conjugate momentum of x. This system is a dissipative system and it is easy to 

derive 

dM o , , 

The equahty in (||) is hold only y — 0. Because y — and a; 7^ is not a stable point, H decreases. 
Since the potential V{x) has a minimum only at a; = 0, Eq. (|9|) gives 

Hiy,x), y(a:)^0 (10) 

in the limit oi t —^ 00. Next, we define 



r = yx^T^ (11) 

In the large t limit, r^ also approaches zero, since 

^r^ < H{y, x) < max Q, ^ + ^ ^'- (12) 

Therefore, x and y approach zero in this limit. Similarly, we can show 

r^^oo, (13) 

when t approaches —00. Since a; = and y — x ^ ±cx) is not a stable point, we consider the case 
x^ — > 00 in the limit t -^ —00. In this case, the nonlinear equion (0) is reduced to a linear equation 

3 

x + 2i+-a; = 0. (14) 

4 

This equation can be solved analytically and we get W 

x{t)^Cie'i+C3e-^, (15) 

where Ci and C3 are some constants. To satisfy the infra-red boundary condition (IRBC) (H), 
C3 = 0. We remark that, due to Eq. (o), vi and v^ are conserved in this limit. On the other hand. 



in the high-energy region in which a; — )■ 0, the nonlinear equation (0) reduces to another Hnear 

equation as 

3 
x + 2x+-x + Xx^0 (16) 
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The general solution in this limit is given by Q| 

x{t) = cie-(i-'^)* + C2e"(i+'^)* (A<l/4) 



x{t) = C3e~* + C4te"* (A = 1/4) 
x{t) = C5e~^ cos {bt + cq) (A > 1/4), 



(17) 



where a = ^^1/4 — A, b = y/X — 1/4 and ci ~ cq are some constants. Therefore, in the high-energy 
region, the solution has oscillating behavior in the case of A > Ac = 1/4 and it can connect the line 
vi^O {y = -x/2) of IRBC and the hne ug = (y = -3x/2) of UVBC with mo = 0. In that case, 
the equation (nl) has a non-trivial solution B{p^) ^ and DSB of chiral symmetry takes place. 
& 0, 0| Therefore, ac = 47rAc/3 — 7r/3 is called a critical coupling. Below, for simplicity, we call 
Ac itself a " critical coupling" . 

The condition A > Ac is necessary but is not sufficient for DSB. If the solution continues to 
oscillate in the ultraviolet region, it may go beyond the line y = — 3a;/2 of UVBC. On the other 
hand, the integral (|l|) may diverge. To converge the integral, a form factor (or a cutoff) which 
switches off the interaction at very large t may be introduced. Due to the form factor, the equation 
(H) becomes ( [l4| ) again and the solution ceases to oscillate. If the solution satisfies UVBC when 
the interaction is switched off by the form factor, it evolutes on the line y — — 3a;/2 after that time. 
This is a sufficient condition for DSB. The aim of this paper is to re-examine the nonlinear effects 
in this scenario and show how the linear approximation works well. To achieve this, we define 



X — r cos 9 and y = r sin t 
We get the differential equations in this polar coordinate. 



dr ( 1 

— = sin 9 [ —2r sin 9 H — r cos ( 

dt \ 4 



d9 , „ . „ cos2 6' 



A- 



r cos^ 



1 + r^ cos-^ 



cos 



— = -l-2cos6'sin6l-^ — A- „ ,^. 

dt 4 1 + r2 cos2 9 

The equation (EOh can be rewritten as 

d9 _ ?/ + 2xy + (3/4 + Xeff{x))x^ _ y^ + 2xy + iJ^{x)x^ 
dt r^ r^ 



(18) 



(19) 



(20) 



(21) 



where Xeff{x) = A/(l -I- a;^)(< A) is an effective coupling and uj'^{x) — Jl^' jx = 3/4 + Xc,ff{x). 
The ^ becomes negative unless a second order equation for y 



y"^ + 2xy + (3/4 + Xeff{x))x'^ = 
has real solutions. This condition gives 



(22) 



Ae//(x)<l/4=Ac. 



(23) 



In this case, the real solutions are given by 




^-x± ^x^il - u;^{x)) = -x±Jx^(-- Xeffix) . (24) 



These are boundaries of a region where ^ > 0. We call this region " typhoon- like" because the 
solution has a lef t-ha.nded rotation in this region. On the boundaries ^ = 0. For convenience, 
we include the boundaries in the typhoon-like region. In the remaining part of x-y plane, ^ < 0. 
We call this region "anti-typhoon-like" because the solution has a right-ha,nded rotation in this 
region. It is interesting that two boundaries lie between the line y = —x/2 of IRBC and the line 
y = —Sx/2 of UVBC on x-y plane. In particular, two boundaries coincide with the line y = ~x/2 
and the line y = —3x/2 when A = 0. In fig. 1, we show the boundaries in the cases of A < Ac 
and A > Ac- The typhoon-like region is continuous at the origin of x-y plane, if A < Ac. It is 
discontinuous at the origin, if A > Ac- 

Our picture of DSB is as following. In the infra-red limit, the solution is given by x = Cie~*' ^. 
As t becomes large, the interaction is switched on and the conservation of vi is violated, i.e., the 
solution deviates from x = Cie~*/^. As x becomes smaller, the effective coupling Xeff{x) becomes 
larger. However, if A < Ac, the effective coupling Xeff{x) can not be larger than Ac. Therefore, 
there is a typhoon-like region just between the line y = —x/2 of IRBC and the line y — ~3x/2 of 
UVBC. Since the solution has a left-handed rotation within the typhoon-like region, the solution is 
prevented from connecting the line y — —x/2 with the line y = —3x/2 in the same quadrant. On 
the other hand, the anti-typhoon-like region prevents the solution connecting the line y — —x/2 
with the line y = —3x/2 in the different quadrant. If A > Ac, at some t, \eff{x) becomes greater 
than Ac. The typhoon-like region disappears. The solution begins to oscillate and can connect 
IRBC and UVBC. The DSB may happen. Therefore, A > Ac is a necessary condition for DSB. 

As is seen above, the nonlinear effects can be included in the effective coupling Xeff{x). The 
typhoon-like region also prevents the solution deviating much from the line y = —x/2 of IRBC. 
Therefore, x — Cie~*" may be a good approximation for the exact solution, until Ae//(a;) becomes 
Ac. Furthermore, if A is not much larger than Ac, x becomes sufficiently small before Ae//(a;) 
becomes Ac, and, therefore, the linear approximation (Hq) works well in the region where Ae//(x) > 
Ac. In this case, the effective coupling Ae//(x) can be approximated as A/(l -I- (Ci)^e~*). This is 
nothing but the linear approximation used by Fomin, Gusynin and Miransky. [pi 

Next, we consider a sufficient condition for DSB. Suppose that A is slightly larger than Ac and 
X does not deviate much from y = x = —x/2 until some to where r becomes sufficiently small. We 
can re-define to as an origin of t. Neglecting the r-dependence, we integrate the equation (|o|) with 
an initial conditions x{0) = — a;(0)/2 and impose a final condition x{tcut) = —^xitcut)/'^ where tcut 
is a cutoff parameter. The result gives the condition 

bU^t^nn-e* (n = l,2...), (25) 

where b is the same as in Eq. (^^ and 9* — arctan jrfp- • This is an approximate version of the 
condition for DSB in |]. 

Our analyses can be generalized in more general equation as 

dV(t x) 

X -\- 27(t, x)x -\- uj'^{t, x)x = x + 2-y{t, x)x H ^^-^ = 0. (26) 



In this case, we get 



dx 



dH „ , ^ 2 dVit,x) 

^ = -27(t,.).^ + 4^ (27) 



-dt = 7^ (^^^ 

As is already seen, if we put 7(t, x) = \ and uj{x, t) = 3/4 + A/(l + (Ci)^e~*) in Eq. (p6|), we get 
the linear approximation used by Fomin, Gusynin and Miransky. p, 0| If we put 7(^,3:) = 1 and 
u}{x^ t) = 3/4 + A9(t) where Q{t) is a step function, it reproduces the results in Eq. (|25|). 
To apply our analyses to the improved ladder calculation, consider the integral 

\{z,u)B{u)u f°° X{z,u)B{u) 



z Jo u + B^u) J^ 



du ,' D2^ ^ +™o, (29) 



where z ^ p^ — m^e* and u ~ q^ . For A(a;, y), we use the Higashijinia-Miransky approximation 

1,1 

A(z, u) = \{z)Q{z -u) + \{u)Q{u - z), (30) 

Differentiating B{z) twice with respect to z, we get the differential equation as Eq. ( pq ) with ||, |l[ 

7(i,x)=7(t) = l-^C/(t); f^(0^^ (31) 



and 



2, , I9y(i,x), ^^, .1/3 l^^.A 2 A*(i),_,, , 2^ 



^'(i,^) = -^^; m2:)^-(---C/(i))a;^ + ^^^log(l + a:^), (32) 



where A*(t) = A(t) - \{t). The IRBC and UVBC are modified as 

(i ^ -c») (33) 



3t (i + ^x) 



^~ \*{t) 

and 

We remark that IRBC ( p3| ) is reduced to the original one Wi = if \*{t) is finite in the limit 
t -^ -00. Due to (|23), if 



D = x\^\t)-u\t,x)) = x^ ( 1(1 - C/(i))^ - Kffit.^)) ; A:^/(i,x) ^ ^1^ (35) 

becomes negative, a typhoon-hke region disappears. The condition ( |35| ) gives X*ff{t,x) > 1/4(1 — 
U{t))'^. In the case of Z? > 0, the boundaries of the typhoon- like region is given by 



y = -xj{t) ± VD. (36) 

If A(t) and X{t) is smah compared with X{t), we get 7(i) « 1, w^(a;, i) « 3/4 + A(i)/(1 + x^) and 
Ug « W3. In this case, approximately, max(A(i)) > Ac is a necessary condition for DSB. 

In summary, we have re-examined the solution of DSB of chiral symmetry in the ladder SD 
equation. Without further approximation, we have showed that the nonlinear effects are included 
in the effective couplings Xeff{x). If Ae//(a;) > Ac, the typhoon- like region disappears and DSB 
takes place. The typhoon-like region prevents the solution deviating much from IRBC and makes 
the linear approximation works well. The analyses is generalized in the case of the improved ladder 
calculation with the Higashijima-Miransky approximation. 
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Figure 1: The typhoon- like region (denoted by "TP") and its boundaries in x-y plane. The 
remaining part in the plane is anti-typhoon-like. Two straight lines are IRBC and UVBC with 
mo = 0. (a) A < Ac. (b) A > Ac. 



